Swimmers in thin films: from swarming to hydrodynamic instabilities 
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We investigate theoretically the collective dynamics of a suspension of low Reynolds number swim- 
mers that are confined to two dimensions by a thin fluid film. Our model swimmer is characterized 
by internal degrees of freedom which locally exert active stresses (force dipoles or quadrupoles) on 
the fluid. We find that hydrodynamic interactions mediated by the film can give rise to spontaneous 
continuous symmetry breaking (swarming), to states with either polar or nematic homogeneous 
order. For dipolar swimmers, the stroke averaged dynamics are enough to determine the leading 
contributions to the collective behaviour. In contrast, for quadrupolar swimmers, our analysis shows 
that detailed features of the internal dynamics play an important role in determining the bulk be- 
haviour. In the broken symmetry phases, we investigate fluctuations of hydrodynamic variables of 
the system and find that these destabilize order. Interestingly, this instability is not generic and 
depends on length-scale. 



In nature, baths of micron-scale swimmers are found to 
show remarkable out of equilibrium phenomena ranging 
from anomalous diffusion and viscosity enhancement, to 
turbulent and swirl-like motion or self organization into 
complex dissipative structures [IHH . They act as the in- 
spiration for man-made devices able to control and mix 
fluids on micron scales. The search for the design prin- 
ciples of such devices remains in its infancy. How- 
ever, simplified models of low Reynolds number swim- 
mers should prove to be a useful starting point for the 
theoretical understanding of this class of collective phe- 
nomena. 

A very simple picture Q then, of a swimmer is an inter- 
nally driven vector oriented towards its direction of mo- 
tion. Theoretically the dynamics of a collection of such 
objects can be described, on long length- and time-scales, 
by vector and tensor equations which are natural gener- 
alizations of equilibrium liquid crystalline hydrodynam- 
ics (To| . These active fluids involve the study of conserved 
and broken symmetry variables that are non-equilibrium 
analogues [Tl| of Goldstone modes. Previous studies 
on active suspensions [12Nl5j| have shown that fluctua- 
tions in these modes, destabilize ordered states. In 3d 
this instability is termed generic, as it is independent 
of length-scale. However, these descriptions are in fun- 
damental sense, vhenomenological since another physical 
must be invoked to generate the or- 



mechanism 14 



, phe 

M 



dered states which have subsequently been shown to be 
unstable 



12|. 



In this letter, we introduce and study a self-contained 
and microscopically defined physical system in which it is 
possible to both generate homogeneously ordered states 
and to examine their stability. This is provided by a 
'suspension' of swimmers confined to two dimensions by 
a viscous thin film In addition the thin film geom- 
etry is particularly accessible to experiments both from 
the point of view of ease of observation and external acti- 
vation. We study analytically a model that directly links 
the collective behaviour to the microscopic dynamics. It 



combines fluctuations, both active and passive, with the 
deterministic motion due to activity and hydrodynamic 
interactions. 

We find that the purely physical coupling mediated by 
the thin film can give rise to local spontaneous breaking 
of symmetry and we identify the possible ordered states. 
The system can have both nematic order ; characterized 
by a macroscopic axis of mean orientation n and symme- 
try n — > — n, or polar order ; with mean orientation axis 
p for which p ^ -p. We then examine the stability of 
each of the homogeneous phases to hydrodynamic fluctu- 
ations. We find that the isotropic phase is stable. On the 
other hand, as for the bulk 3d system, the homogeneous 
broken symmetry states are destabilized by the hydro- 
dynamic modes of the system. Here, however the thin 
film weakens this effect and the instabilities observed are 
length scale dependent (i.e. not generic). 

We restrict ourselves to a dilute solution, and we 
make use of a mean field approximation and study 
the one-particle distribution function: c(R, u, t) = 
— Tl)5(ui — u)}, the probability of finding a 
swimmer with average orientation u at position R. This 
satisfies a dynamic equation 

d t c = -V -St-H-Jr ; (1) 
J T = -D T • Vc+ (vu + V)c, J R = -D R TLc + Sic, 

in terms of translational J<r(R, u, t) and rotational 
^Tr(R, u,t) currents, with U := ti A ^ 0. Dy and 
Dr are respectively the translational (rotational) diffu- 
sion tensor (constant) of the swimmer and represent both 
passive and active fluctuations. The deterministic quan- 
tities are active, describing self-propulsion: vu, and the 
translational V and rotational f2 velocities induced on a 
swimmer due to the activity of others by hydrodynamic 
interactions mediated by the film. 

Thin film hydrodynamics: The film is described as an 
infinite incompressible two dimensional layer of fluid with 
(2d) viscosity n filling the plane z = and coupled hydro- 
dynamically to another incompressible bulk fluid of (3d) 
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viscosity rj e which fills the region z ^ 0. To distinguish 
between them we indicate the three-dimensional quanti- 
ties with a prime. We consider the fluid dynamics in the 
vanishing Reynolds number (Stokes) limit where inertia 
can be neglected [l?}- For the in-plane quantities, given 
an in-plane force densit y F {x, y) the velocity v(x, y) and 
pressure p(x, y) satisfy [18| 



r]V 2 ± v + W±p + cr+ - ol 



F: V, 







(2) 



where Vj_ = (d x ,d y ) is the 2d gradient operator, af := 
n e d z v'\ ± is the shear stress of the bulk fluid at the 
top/bottom of the thin film; in the external region, z =/= 0, 
the velocity v'(x, y, z) and pressure p'(x, y, z) satisfy the 
Stokes equation : ?7 e VV + Vp' = ; V • v' = , 
where V = (V±,d z ) is a 3d gradient operator. The ratio 
of the two and three dimensional viscosities introduces 
a length-scale s := ij/{2rj e ) that governs two asymptotic 
regimes. For r < s dissipation occurs almost entirely in 
plane and hydrodynamic flow fields are quasi-two dimen- 
sional while at lengths r»s dissipation is mostly due to 
flow out of plane, and the hydrodynamics is similar (but 
not identical) to that in three dimensions 

The viscous drag coefficient on a flat disk of radius a 
embedded in the film is 7 = Airrj/g, with g a function 
of s/a [19 ]. On length scales large compared to a, the 
interaction between several disks lying in the film can be 
approximated using point-like forces at their centers and 
the Green's function H of eq. ©, corresponding to the 
flow, v(r) = H(r — ro) • fo generated by a point-like force 
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fo at Tq. The tensor H(r) 
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the thin film equivalent of the Oseen tensor 
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In the 

following, we work in the limit r s 

We consider swimmers with an average speed v mov- 
ing in a direction u which can be represented at large 
length-scales as a time-dependent force dipole that gener- 
ates an associated velocity field. We also consider swim- 
mers for which the force dipole is zero and the leading 
behaviour is determined by a force quadrupole. Hence 
the average force density of a swimmer with mean posi- 
tion, r a is of the form f Q (r) = — f d LuiiiV i5(r — r Q ) + 

\f q L uuiUjVi\7jS(r — r a ) H , representing the dipole 

and quadrupole contributions respectively where L is the 
typical dimension of the swimmer. 

For a concrete calculation, an explicit microscopic 
model of a swimmer is required and we have used a 
three-disk model of swimmer [2l| . This is character- 
ized by a minimum number (2) of degrees of freedom, 
with typical length I that move in a non-reciprocal fash- 
ion in time with frequency oj = to achieve loco- 
motion in the Stokes limit. In the following we indi- 
cate the time-average over a swimmer cycle period T 
with an over-bar h = ij; h(t)dt. In the limit of 
small sinusoidal oscillations of amplitude d around I is 
possible to obtain [2l[ the average self-propulsion ve- 



locity as 
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r] where a := |f 



and our convention is that the swimmer move in the 
same direction of it. For the average force we find 



L = I and f d 
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/ := 7wd is the force scale on each degree of freedom 
of the swimmer. The parameter k, which is the average 
ratio of the internal lengths, controls the nature of the 
swimmer: the dipole is positive (pusher) for k > 1, neg- 
ative (puller) for k < 1 and zero (quadrupole) for n = 1. 

The interaction between two such swimmers, which is 
the origin of V, fl in cq. (fT|), is complex as each one 
is characterized by periodic internal dynamics whose cy- 
cles may also have different phases [22j. Here we will 
restrict ourselves to swimmers that all have the same 



phase [22|, and the limit where the typical separation 
r between the swimmers centres R A and R B is much 
larger than the typical dimension I of the object. Dy- 
namical quantities depending on hydrodynamic interac- 
tions, such as translational V26 := R Q — vu a and an- 
gular u)2b '■= u a A u a velocities felt by each swimmer, 
for a — A,B, are described to good approximation by 
the first few terms of the expansion in spherical harmon- 
ics [l5, 22 1 of the tensor H. To leading order, for swim- 
mer A we find v 2fc = f^{(3£ 2 - 1)£ -B^} + 0(l/r 3 ) 
se<2> {3(^ + 2BC - 5AB 2 ' -' 

4\ tT„„„ A — (A . S.A 
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and u)2b 

C)H^vuf } + {l/r A ). Here A:=(f- u A ), B := (r 
C := (u A ■ u B ) where u A and u B are swimmers orienta- 
tions, r is the separation vector and is related to the 
time-dependent force dipole and scales as fl. We neglect 
the effect of interactions on the internal dynamics and so 
do not address synchronization effects [20(. 

The collective dynamics is obtained by a coarse grain- 
ing procedure, first in time and then in space. The ve- 
locities f2, V in eqs. ([1]) are obtained from the two body 
velocities for swimmers with positions R and R' and ori- 
entations u and u' as 
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(R — R', u, m')c(R', u') 
(3) 



where, as before, the over-bar denotes time average over 
the swimmer cycle. Order parameters such as local den- 
sity p, polarization P and ncmatic orientation tensor S 
are defined as moments of c : 



(R,t) 




1/2), 



c(R,u,t). (4) 



Homogeneous ordered states These are states in 
which c and its moments defined above do not vary 
with position (denoted by c° , p° , P° , S° respectively). 
Under these conditions, the mean field velocities are 
n° = /3 u A P° + Pxu A S° ■ u and V° = f3 T P° where 
the coefficients /3q, Pi and fix are averaged quantities 
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that dep end on the microscopic details of each swim- 
mer |23l |. To leading order /3 « / ngpA 3 fj j( k2 ~ -*-) an< ^ 

They are both positive for k > 1 (pushers) and negative 
for k < 1 (pullers). Inserting the expressions of $1° and 
V° in eq ([T]) , taking the time derivative of eq Q , we 
obtain dynamic equations for the moments of c . Density 
is conserved, hence dtp = 0. The others are, 



As in many active systems, we find that their effect 
is to destabilize order 



d t P° 
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•°; (5) 



Analysis of eqs © shows that, when /3 or /?! are pos- 
itive, the system can undergo a bifurcation, that sig- 
nals the appearance of order. Hence pushers can give 
rise to order whereas pullers cannot. The conditions 
-D R + p°A)/ 2 = and -AD R + p°/3 1 /2 = define two 
critical lines in the space of parameters p° and /, above 
which the instability occurs 24|. The ratio /3i/(4/3q) 
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FIG. 1: Phase diagram showing homogeneous states and their 
hydrodynamic stability. Isotropic-Nematic and Nematic- 
Polar transitions for dipolar swimmers. The transitions can 
occur only for pushers. Force / is measured in arbitrary units. 
Parameters are set to - = Dr = 1, d = 3ao, / = 3d, A = 2, 
and k = 1.2. In the inset is shown the I-N transition for 
quadrupolar swimmers (same parameters except « = 1) . 

determines if the I-N or I-P transition occurs at lower 
density (see Fig. [1] where we have chosen parameters 
for which I-N occurs first). For pullers (3q and Pi are 
negative and contribute to enhance the noise (diffusion) 
in the system. For quadrupolar swimmers, these lead- 
ing order terms vanish and in addition to higher order 
terms one must keep track of internal mode dynamics of 
the swimmers. The leading order terms in /3o are nega- 
tive, so enhance diffusion while in j3\ they are positive 
and promote order. Hence, for this particular micro- 
scopic prescription, we find no polar phase. We find 

ft « ffj nfmxsi 1 + f( 27r + !)] & A P lot of thc 

critical line in this case is shown in the inset of Fig. [T] 
Clearly, this transition occurs at higher densities than 
dipolar swimmers since it is due to higher order terms. 

We now discuss fluctuations in hydrodynamic variables 
about isotropic and ordered states for dipole-swimmers. 



121. 



In the following we con- 
sider the deviations of the fields from their homogeneous 
values given by Sp = p — p°, similarly for P and S. 
We introduce Fourier transforms in the usual way as 



/(k),/(k)= r r e ik -/(r). 



/(') = Jk 

Isotropic state In the homogeneous isotropic state 
p = p° , P° = S° = and p is thc only hydrodynamic 
variable. Variables SPn := SP ■ fc, splay 8S\\\\ and bend 

(i5S||_l) q := 5Sbck c {Sab — k a kb) show diffusive behaviour. 
From the resulting set of equations density fluctuations 
are linearly stable [TJ] ■ 

Polar state In the homogeneous polar state p — p , 
P = P° and the hydrodynamic variables are p and the 
director P. The orientation tensor is slaved to P, and 
given by S = SpP 2 (PP — |) where Sp is determined 
by eqn. ([5]). The magnitude P is not a hydrodynamic 
variable and relaxes to a constant value on microscopic 
time-scales. We set P = 1 in the following and study 
linear perturbations around this state. Decomposing 
k = cos4>P° + sm(j)k±, with k± ■ P° = 0, to leading or- 
der in k we find that fluctuations in density and director 

are coupled via splay (k± ■ SP) , giving a growing mode 

T with real part U(T) - fc s5 ^ (2> cos 20(2 + cos 20) and 
imaginary part, which determines thc propagation speed, 
9f(r) - ±k^=\ sin 0|, as shown in figure^a). When 9?(r) 
is positive fluctuations grow exponentially, signalling an 
instability of the ordered state due to hyd rody namic in- 
teractions, analogous to those found in 3d [12j. However 
unlike those, the growth rate here scales as k in the k — > 
limit. 

Nematic state A similar analysis can be performed 
around the homogeneous nematic state (finite p°, §°). 
In this case hydrodynamic variables are density, p and 
nematic orientation tensor, S = S(n ®h — |). Again the 
magnitude S relaxes fast and in the following will be set 
to S = 1. To leading order in k fluctuations of density Sp 
and nematic director Sn are decoupled and splay fluctu- 
ations have a real mode T ~ fc cos 20(2 + p° cos 20). 
For angles above 7r/4 destabilize the order, as shown in 
Fig-IH(b). For quadrupolar swimmers the relevant terms 
describing hydrodynamic interactions scale as fc 2 instead 
of k. Their analysis is not reported here (2p| . 

It is instructive to compare our results with previous 
studies of swimmers in a 3d fluid |12l |l4j, |15( . Our analy- 
sis starts from a microscopic model, and 'integrates out' 
the fluid degrees of freedom to see the effect on the other 
hydrodynamic modes. Alternatively one may perform 
a phenomenological analysis of an ordered state as in 
[l2| replacing the Stokes equation by eq ([2]); this yields 
qualitatively the same hydrodynamic instabilities that we 
have presented above. In essence, the thin film model 
changes the Fourier spectrum of the hydrodynamic ker- 
nel from k 2 to fc + sfc 2 which in the limit sfc 1 reduces 
to « fc. This is the origin of the different scaling of the in- 



4 




research studentship. TBL acknowledges the support of 
the EPSRC under grant EP/G026440/1. 



2 4 6 8 10 

* (rad) 

(a) 



FIG. 2: Angular dependence of the growing mode V in the 
splay fluctuations around homogeneous order for pushers. A 
positive sign indicates instability, (a) polar state. 5ft(F) is 

measured in units of k(sSp\Q \)/(8>v) an d 9(r) in units 
of kv/y/2; (b) nematic state. 3?(r) is measured in units of 
fc(s|9 (2) |)/(4j7) and we have used p° = 0.3. 



stability here. In another contrast, hydrodynamic inter- 
actions of simplified far-field models of swimmers cannot 
lead to homogeneous order (swarming) in 3d 15[. This 
can be ascribed to a mathematical cancellation that oc- 
curs performing angular integrals of spherical harmonics 
generated by the hydrodynamic kernels. In the thin film 
limit considered here, this is circumvented by the con- 
finement of the swimmer directors to two dimensions. In 
addition the instability generated by activity in the film 
is 'soft' in the sense that it scales with k — > in compar- 
ison to the 3d case where it is independent of k (hard) to 



leading order [12j. We can conclude that in comparison 



to 3d, hydrodynamic interactions between swimmers in 
a thin film favour order but are still not strong enough 
to overcome the destabilizing effect of activity. Finally 
we note that our approach is complementary to models 
where simple phenomenological rules of interaction be- 
tween swimmers are used to study aspects of collective 
behaviour [ll, 24, 25 1. In contrast, here we 'derive' such 
rules from a particular microscopic model using a coarse- 
graining procedure which is valid under a precise set of 
conditions. 

The question of how to characterise the system be- 
yond the instability ljj, [25J remains open. This high- 



lights a fundamental theoretical issue. The instabilities 
of homogeneous ordered states seen here and in other 
non-equilibrium active systems illustrate the limits of the 
Landau- Ginzburg framework [l2, 15 1, which has been so 
successful in the study of phase transitions in equilib- 
rium systems. It might be that alternative approaches, 
such as describing the homogeneous state as dynamical 
rather than stationary, or the development of a formalism 
in which hydrodynamic and non-hydrodynamic variables 
are treated on the same level, have to be considered. 
ML acknowledges the support of University of Bristol 
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